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Numerical studies implementing different versions of the collisionless Ohm's law have shown a reconnection 
rate insensitive to the nature of the non-ideal mechanism occurring at the X line, as soon as the Hall effect 
is operating. Consequently, the dissipation mechanism occurring in the vicinity of the the reconnection site 
in collisionless systems is usually thought not to have a dynamical role beyond the violation of the frozen-in 
condition. The interpretation of recent studies have however led to the opposite conclusion that the electron 
scale dissipative processes play an important dynamical role in preventing an elongation of the electron layer 
from throttling the reconnection rate. This work re- visits this topic with a new approach. Instead of focusing 
on the extensively studied symmetric configuration, we aim to investigate whether the macroscopic properties 
of collisionless reconnection are affected by the dissipation physics in asymmetric configurations, for which the 
effect of the Hall physics is substantially modified. Because it includes all the physical scales a priori important 
for collisionless reconnection (Hall and ion kinetic physics) and also because it allows one to change the nature 
of the non-ideal electron scale physics, we use a (two dimensional) hybrid model. The effects of numerical, 
resistive and hyper-resistive dissipation are studied. In a first part we perform simulations of symmetric 
reconnection with different non-ideal electron physics. We show that the model captures the already known 
properties of collisionless reconnection. In a second part, we focus on an asymmetric configuration where the 
magnetic field strength and the density are both asymmetric. Our results show that contrary to symmetric 
reconnection, the asymmetric model evolution strongly depends on the nature of the mechanism which breaks 
the field line connectivity. The dissipation occurring at the X line plays an important role in preventing the 
electron current layer from elongating and forming plasmoids. 



I. INTRODUCTION 

Magnetic reconnection is a universal plasma phe- 
nomenon releasing magnetic energy into thermal and 
bulk kinetic energyi~— . Since reconnection is thought 
to be at the root of very fast energy release events, a 
major issue is the understanding of the physical mecha- 
nisms controling the rate at which the magnetic flux is 
being reconnected, and at which scale they operate. Re- 
connecting field lines requires that all the plasma species 
violate the frozen-in condition. In collisionless systems, 
this leads the structures to thin down to scales where 
finite Larmor radius effects become important. The elec- 
trons are the lightest species, they will therefore control 
the scale at which reconnection ultimately occurs in such 
systems, but not necessarily the scale that controls the 
overall rate. The puzzle then consists in understanding 
whether the electron scale physics, the (larger) ion scale 
physics or a coupling between the two, controls the over- 
all dynamics of the reconnection process and how. 

Modeling magnetic reconnection with different simpli- 
fied versions of the collisionless Ohm's law has led to 
great progress in understanding which mechanism plays 
a pivotal role in the collisionless reconnection process. It 
is now generally accepted that the Hall physics plays a 
crucial role in enabling a fast regime, i.e. a rate that 
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is large enough for reconnection to be consistent with 
the phenomena it is considered to be a key ingredient 
of. Multiple numerical studies have indeed shown that 
including the Hall effect in reconnection models results 
in a large rate, which is moreover quite insensitive to the 
non-ideal electron dynamics occurring close to the recon- 
nection site^'^. 

A paradigm, generally called the Hall reconnection 
models has been proposed to explain why the Hall effect 
alone results in these two facts^i^— . First, the Hall ef- 
fect, owing to its dispersive properties, is thought to make 
the outward electron flux insensitive to the mechanism 
breaking the field line connectivity. A Sweet-Parker— 
like reasoning then leads to the conclusion that the re- 
connection rate must remain unaffected by any parame- 
ter altering the physics at this scale, no matter what it 
actually is as long as it falls well below the ion inertial 
length. The second consequence of the Hall effect is the 
propagation of Hall fields along the separatrices, which 
opens the exhaust by heating the ions and accelerating 
them to maintain a large downstream flux, independent 
of the system size. 

Understanding whether collisionless reconnection is 
sensitive or not to the electron scale physics occurring 
at the X line or is entirely controlled by the Hall ef- 
fect alone, however remains an open question. Recent 
simulations ^^'^^ of domains larger than before and/or 
with open boundaries, have revealed an elongation of the 
electron current layer in the downstream direction. Ev- 
idence was found to support the idea of a relationship 
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between this elongation and the concomitant decrease of 
the reconnection rate. Two mechanisms were proposed 
to explain the cessation of elongation. One mechanism is 
the triggering of a secondary tearing instability, which, by 
breaking the current sheet into multiple plasmoids, would 
dynamically prevent further elongation of the layer, and 
at the same time, make reconnection a highly unsteady 
process. A second mechanism would be the balance in 
the downstream direction of the outward propagation of 
magnetic flux with the dissipation occurring at the X 
line, which could occur for long periods of time and make 
the reconnection process appearing as quasi-steady. The 
question of what would cause an unbalance between these 
two processes and then result in the elongation of the cur- 
rent sheet has however not been addressed. 

On the other hand, other fully kinetic simulations^iiS, 
in quite similar initial conflgurations, have shown a fast 
and steady reconnection rate and insensitive to variations 
of the electron mass, supporting the Hall reconnection 
paradigm. As a result of these simulations, it was con- 
cluded that the extent of the electron layer in the down- 
stream direction is controlled by the Hall effect alone 
and its ability to rotate the reconnected magnetic field 
in the out-of-plane direction with a speed depending on 
the scale. Furthermore, Hall MHD simulations without 
explicit dissipation physics, i.e. with only numerical dis- 
sipation, have revealed a fast reconnection rate, similar to 
the one measured in fully kinetic systems, and moreover 
independent on the mesh resolutioni^ii^. These last re- 
sults are consistent with the electron dissipative physics 
playing no role and give the Hall effect the full control of 
the reconnection dynamics. 

Finally, recent theoretical investigations^^ have shed 
new lights on the physical role of the electron dissipa- 
tion scale mechanisms. The structure of the electron full 
pressure tensor at the X line has been explained from 
the kinetic viewpoint to be related either to the electron 
bounce motion around reversed field lines, or to the mix- 
ing of cold/slow electrons entering the layer with hot /fast 
electrons leaving it, depending on the coplanarity of the 
macroscopic magnetic configuration. The theory showed 
that this non-gyrotropic contribution to the electron mo- 
mentum can be related to a viscous effect from the fiuid 
viewpoint, the off-diagonal components of the pressure 
tensor being linked to the second derivative of the elec- 
tron bulk velocity. Physically, the role of the reconnec- 
tion electric field can thus be understood as the way to 
sustain the current density required at the X line to main- 
tain the magnetic shear, which would otherwise be dis- 
sipated by the diffusion of electrons accelerated in the 
layer and leaving it. 

In this paper, we aim to study further the sensitivity 
of collisionless reconnection to the dissipation physics. 
Since the classic reconnection setup seems to lead to 
rather different conclusions, we adopt a different ap- 
proach, and decide to investigate how sensitive collision- 
less reconnection is to the dissipation physics in a more 
general, asymmetric, configuration too, i.e. where the 



reconnecting current sheet is no longer separating two 
identical plasmas, in terms of density and magnetic field 
strength. The initial perfect symmetry of the current 
sheet, extensively used to initialize numerical models, is 
indeed hardly expected anywhere, but, to some extent, 
at the Earth magnetotail. Besides, the few studies which 
have been so far focused on asymmetric reconnection 
have shown that the structure of the Hall fields, on which 
relies the Hall paradigm, is significantly modified^ —. It 
is therefore important to understand how these systems 
respond with regard to the current debate on the role 
of the dissipative physics. We will focus on the role, 
in reconnection, of the dissipation, as defined from the 
fiuid viewpoint, as any mechanism that is not consid- 
ered as ideal and which enables field lines to change their 
connectivity. Previous understanding may lead to think 
that any dissipation effect coming after the Hall term 
in the Ohm's law, should leave the overall evolution of 
the reconnection process unchanged. This question does 
not need a fully kinetic code to be addressed, as the hy- 
brid kinetic formalism, includes all the ingredients previ- 
ously mentioned as part of the steady Hall reconnection 
paradigm, neglects electron non-ideal effects to which the 
paradigm attributes no dynamical role and let us control 
the dissipative physics by choosing different mechanisms 
and investigate their respective role. Any deviation from 
the standard steadiness and insensitivity to the dissipa- 
tion mechanism would emphasize the need for a better 
understanding of reconnection and would allow one to 
identify key ingredients for its modeling. Furthermore, it 
is important to understand to what extent the use of a 
dissipation mechanism, versus another one, in a hybrid 
code, changes the macroscopic evolution of asymmetric 
reconnection, as these codes are the only one able to sim- 
ulate large scale phenomena and include the appropriate 
ion dynamics at the same time. This paper does however 
not investigate whether the chosen dissipation models 
can be justified by specific underlying kinetic processes, 
and only looks, from the macroscopic viewpoint, whether 
any model results in the same overall evolution of recon- 
nection or not, in symmetric and asymmetric configu- 
rations. The comparison of that overall evolution with 
a fully kinetic model is a topic explored in a separate 
paper—. The investigation of how the kinetic processes, 
occurring at the X line in asymmetric reconnection, can 
result, from the fiuid viewpoint, in dissipation, and how 
it compares to the present dissipation mechanisms is also 
considered off-topic and should be investigated in the fu- 
ture using, this time, fully kinetic simulations. 

In the second part of this paper, we will present in 
detail the numerical model used to perform this study. 
Then, the third part will be focused on hybrid simula- 
tions of symmetric reconnection with different dissipa- 
tion mechanisms. This part will show to what extent our 
model recovers the already known properties of Hall re- 
connection. In a fourth part, we will investigate the role 
of dissipation in asymmetric reconnection. Finally the 
fifth part will summarize and discuss our findings. 
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II. NUMERICAL MODEL 

The hybrid model we use solves the ion (assumed to 
be protons) kinetic dynamics using the particle in cell 
method^i. Their motion is obtained from equation ([1]), 
where is the proton mass, E and B are the electric 
and magnetic field, respectively. This equation is solved 
using the Boris algorithm^^. The electron population is 
assumed to be a fiuid, meaning that the infinite chain 
of moments of their distribution function is truncated at 
the pressure Pg order, for which an isothermal law is as- 
sumed ([5j) where Tg is the constant and uniform electron 
temperature. The electron bulk inertia is also neglected. 
At all locations and times, the electron fiuid is assumed 
to have the same density as the proton density n, so 
that quasineutrality holds. The proton density itself is 
obtained from the proton velocity distribution (|6|). The 
electron momentum equation is used to calculate the elec- 
tric field in the form of an Ohm's law ([3]), where R is a 
dissipation term. Consistently, the displacement current 
is neglected in Maxwell- Ampere's equation (|2]), where the 
definition of the current density j = en (v^ — Vg) has been 
used, with e the elementary charge, and Vg are the 
proton and electron bulk fiows, respectively. Similarly to 
their density, the proton bulk fiow is calculated from the 
velocity distribution ([7]). The magnetic field is evolved 
using Faraday's law (j4]). 



m,^ = e(v^, xB + E) (1) 

V X B = /ioen (v^ - Vg) (2) 

E = -v^ X B + — (j X B - VPg) + R (3) 

ne 

- = -VxE (4) 

Pg = nksTe (5) 

^ = y f {T^,^pi)dvpi (6) 

Vi = - / ^rpif (r, Vpi) dvpi (7) 



The dissipation (term R in eq. ([3])) is alternatively 
chosen to be negligible (R = r^j ^ 0) in a way that it 
is not able to prevent the collapse of the current sheet 
thickness down to the grid scale but sufficient to pre- 
vent the simulation from crashing due to accumulation 
of noise, resistive (R = T^j) and non- negligible so that the 
current sheet thickness is controlled by Joule diffusion, or 
hyper-resistive (R = — i^V^j). The resistivity r] and the 
hyper-resistivity u are both considered constant in time 
and spatially uniform. Classical resistivity is related to 
electron-ion collisions, whereas hyper-resistivity can be 
seen as a model of electron viscosit}^. Other models 
for non-ideal terms exist but we have not implemented 
them for the following reasons. Localized resistivity is 
known to be the source of fast reconnection without the 
Hall tern>2^ and thus seems inadequate for this study. 



Moreover, unless a somewhat arbitrary and complicated 
model is chosen, it has a constant and not self-consistent 
scale and is spatially fixed, which can be problematic 
considering the possible drifts of the X line seen in asym- 
metric system^ii^i^. Electron bulk inertia can also be 
non-negligible around the reconnection site. Usually, this 
term is implemented for numerical reasons in two-species 
models because it conveniently changes the whistler dis- 
persion relation and prevents whistler waves to severely 
limit the CFL condition. However its standard imple- 
mentation requires one to assume a uniform electron iner- 
tial length where it should be locafii^. Accordingly, it is 
unclear whether the electron mass appearing in such im- 
plementation has the same role as in real systems where it 
moreover also plays a role in the full pressure tensor. Let 
us note, on the other hand, that hyper-resistivity, being a 
viscous operator, shares remarkable similarities with the 
current understanding of the dissipative nature of sym- 
metric electron current layer—. It therefore appears as a 
reasonable and simple phenomenological model, although 
a more general, asymmetric, theory is needed. 

Numerically, the equations for the electromagnetic 
fields are discretized with second order finite differences. 
The magnetic field components being defined at the same 
location, shifted by half a cell from all electric field 
components. Their time evolution is calculated with a 
predictor-corrector scheme. The mesh size is constant 
and uniform, and chosen so that the smallest physical 
scale in ([3]) is resolved. The time step is chosen so that 
the proton cyclotron period and the whistler propagation 
speed are resolved. Considering the thermal and bulk 
speeds in the presented simulations, these criteria guar- 
antee, in particular, that no particle crosses an entire cell 
within a single time step. The results are presented in the 
plane (x, y) and in dimensionless units. The particle den- 
sity and the magnetic field are normalized by arbitrary 
quantities no and Bq respectively. Distances are normal- 
ized by the proton inertial length Si = Va/^ci^ where Va 
and Qci are the Alfven speed and the cyclotron frequency 
based on the density no and magnetic field Po- The re- 
connection of magnetic field is initialized with a localized 
magnetic perturbation. For all the runs presented in this 
paper, the boundaries are periodic in the x direction, and 
closed and perfectly conducting in the y direction. 



III. SYMMETRIC RECONNECTION 

In this section, our goal is to show how symmetric 
reconnection is affected by the nature of the dissipation 
mechanism and to what extent these findings agree with 
previous ones before studying the more complicated case 
of asymmetric reconnection. To achieve this goal, we 
initialize our code with a symmetric tangential current 
layer. The magnetic field has initially only one compo- 
nent Bx {y) = tanh ((^ — ^o) / A), where yo is the middle 
of the domain in the y direction, and A = 0.5 is the half 
thickness of the magnetic field reversal. The normalized 
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density is chosen to be uniform and its value is set to 
unity. The asymptotic ratio of the thermal pressure and 
magnetic pressure /3 is 1. The proton temperature is 
chosen to be isotropic everywhere and is obtained from 
the pressure balance condition nT + = est. The 

electron temperature is set to Tg = 0.25, is constant in 
time and spatially uniform. The protons are initially 
loaded as locally Maxwellian velocity distribution func- 
tions. This initial condition is not a Vlasov-Maxwell 
equilibrium and small magnetosonic waves are emitted 
from the current sheet within the first proton cyclotron 
period. These waves might modulate a bit the reconnec- 
tion rate as they are reflected back to the current sheet 
from the closed top boundaries, but do not contribute to 
the presented results. Unless mentioned otherwise, the 
domain size in the x direction is Xm = 250 and ym = 50 
in the y direction, so that recirculation of plasma and 
waves do not affect the results within the time of interest. 

We will study the results of four simulations differ- 
ing only by the dissipative physics. The first two runs, 
labelled S^^ and S^^^ aim to be cases where the Joule 
dissipation is not strong enough to prevent the current 
sheet thickness to collapse down to the grid scale. The 
number of cells of run S^^ in each direction is set to 
{fix^riy) = (2500,500). The spatial resolution is then 
(Ax, Ay) = (0.1, 0.1). The run S^^ has the same number 
of cells but the size of the domain in both directions is 
doubled, so that the mesh resolution is 0.2 in both di- 
rections. This will allow us to see whether the resolution 
has an impact on the result in these dissipationless cases. 
We have used 7^ = 2 10"^. As a result of such a small 
resistivity combined with these mesh sizes, the current 
sheet collapses at the grid scale, meaning that doubling 
Ay results in dividing the current density by a factor of 
2 (the peak current density being given by jz ~ 1/Ay). 
In this discrete simulated system, the breaking of field 
line connectivity is ensured by finite mesh effects and the 
thickness of the current sheet is obviously limited by the 
mesh size. 

In the third run Sr^^ all parameters are the same but 
the resistivity has been increased up to = 1.4 10~^. In 
the fourth run, we neglect the classical resistivity and 
choose the dissipative electric field to be R = — z^V^j, 
with u = 3 10~^. In this last run, the resolution is 
doubled by increasing the number of cells by a factor 2 
in the y direction. The time step is set to At = 0.001 for 
all runs except run S^^ for which it is At = 0.01. Table 
[J summarizes the run parameters. 

Figures [1] to m show snapshots from the three symmet- 
ric simulations, 6*^^, Srj and at the time they have 
reconnected the same amount of magnetic flux. Run S^^ 
is not shown to save space, but is found very similar to 
S^^ and will be discussed later. On all these figures, the 
top panel represents the run S!^^ {t = 50), the middle 
panel the run Srj it = 57.5) and the run Sjy {t = 52.5) 
is shown on the bottom panel. At first glance the three 



TABLE I. Summary of the parameters 



Name 


V 


V 


Ay 


oa 


2 IQ-^ 





0.1 


qb 

^nd 


2 10"^ 





0.2 


Srj 


1.4 10"^ 





0.1 




2 10"^ 


3 10"^ 


0.05 



runs look very similar, suggesting that the dissipative 
term in eq. ([3]) does not play a critical role in the evolu- 
tion of the system. A closer inspection reveals interest- 
ing details. Fig. [1] shows the out-of-plane magnetic field 
component. As expected for symmetric collisionless re- 
connection, this component has a quadrupolar structure 
around the X point which does not extend upstream of 
the separatrices. This is due to the Hall rotation of the 
upstream in-plane field lines to the out-of-plane direction 
as they cross the separatrices. Looking at cuts in the y 
direction on the top panel of Fig. [TJ one can see that 
if the runs and Sy have very similar structure and 
values (~ 0.4 as in previous works including fully kinetic 
PIC simulationsiiiSi^"— ) , the values of for the run 
Sr^ are somewhat smaller (^ 0.25). 

Figure [2] represents the out-of-plane electron bulk ve- 
locity Ve.z' In the three simulations, v^^z is strongly en- 
hanced at the X point and, to a smaller extent, in the 
separatrix regions. These characteristics are also consis- 
tent with previous studiesi^i^ — . In both runs Snd and 
6*^, the current sheet thickness is comparable to the grid 
spacing in the y direction, whereas it is thicker in run 
Sy. This indicates that, for this resolution at least, the 
resistivity is not large enough to efficiently dissipate the 
incoming magnetic flux upstream of the X point. The 
Joule diffusion seems however to play a role in control- 
ling the length of the layer, as it can be noticed that in 
run S^^ it is close to the x grid spacing whereas it ap- 
pears longer in run Sr^. This last point is more evident 
when observing the in-plane electric field. 

Figures [3] and m show the and Ey components of the 
electric field, respectively. They are zero in the regions 
upstream of the separatrices and are strongly peaked at 
their location. Peak values are 0.1 and ^ 0.5 for E^ 
and Ey , respectively, which is consistent with what is ob- 
served in other studiesi^i^. Both E^ and Ey are weaker 
in the run Sr^ than in the two other runs. In run 6*^^, E^ 
goes up to the X point since the dissipation region has a 
length comparable to the grid spacing, in the other runs 
it starts a bit further and is yet another evidence that dis- 
sipation increases the length of the current layer. This 
feature can also be observed in fully kinetic simulation^. 
In contrast, Ey is not zero in the current sheet except at 
symmetry line itself, and has a strong bipolar variation 
in the y direction. This structure appears in run Sy^ to 
a smaller extent in run Sj^^ but not in run 6*^^, again be- 
cause of the collapse of the electron current sheet down 
to the grid scale in both directions. The small waves seen 
in Fig. [2] for run S^^ can also be observed in both E^ 
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FIG. 1. Top panel : Out-of-plane component of the mag- 
netic field measured as a function oiy ai x — 117 in S^^ (black 
curve), Sn (red curve) and Si, (blue curve). Then, from top 
to bottom : Out-of-plane component of the magnetic field 
{Bz) for run Snd (first panel) at t = 50 , run Sr^ (second 
panel) at t = 57.5 and run Si, (last panel) at t = 52.5. The 
snapshot are averaged over around the indicated time, 
which is when all runs have reconnected the same amount of 
upstream magnetic flux. On each color plot, the black verti- 
cal dashed line shows the position of the cut represented on 
the top panel. 



and Ey components. 

Although the Joule diffusion appears to be not strong 
enough to efficiently dissipate the incoming magnetic flux 
inside the electron current layer, the somewhat smaller 
values of the electromagnetic fields and electron flows 
observed in run Sr^ indicate that it dissipates the energy 
at larger scales. The large scale diffusion is not seen in 
the dissipationless run 6'^^, and not seen either in run Sj,. 
This behavior is consistent with previous understanding 
of the effect of the Joule diffusion in collisionless magnetic 
reconnection^. In run S^^^ small waves can be observed 
upstream of the separatrices. Their absence in the two 
other runs and in fully PIC simulations suggests they are 
a spurious consequence of the lack of dissipation in the 
model. 

Figure [5] shows the reconnection rates measured in all 
simulations. The rate is measured as the out-of-plane 
electric field Ez^ averaged over a unit squared area 
centered on the magnetic flux saddle point, detected at 
each time. The black curve represents the reconnection 
rate of run S^^. Its steady value is End ^ 0.065, 
and goes to ~ 0.13 if normalized to the VAupBup, the 
characteristic electric field at the edge of the averaging 
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FIG. 2. Out-of-plane component of the electron bulk flow 
(vez) for run Snd (top panel), run Sr] (middle panel) and 
run Su (bottom panel). The snapshot are averaged over 
around the same time as in fig[T] 



area. The green curve is the reconnection rate mea- 
sured in run S^^. Appart from the somewhat larger 
fluctuations and the different sampling rate, it is very 
similar and therefore suggests that the mesh resolution 
does not impact the reconnection process dramatically. 
The red dashed curve represents the reconnection rate 
measured in run Su- Its value is Ej, ^ 0.06, slightly 
below End- The rate measured in run Srj (blue line) is 
however significantly below, with a value Er^ ^ 0.055. 
Note that the oscillatory period t e [45, 55] corresponds 
approximately to the time for the waves emitted by the 
current sheet as a consequence of the non-equilibrium 
kinetic initial state, to return back at the center of the 
domain after having been reflected at the top boundaries. 



IV. ASYMMETRIC RECONNECTION 

In this section we change our initial condition to in- 
clude an asymmetry between both sides of the initial tan- 
gential current layer. We choose an initial condition that 
has already been used in previous works^^^^^ for the sake 
of comparison. It consists in a one dimensional tangen- 
tial asymmetric current sheet, where both the magnetic 
field and the particle density balance the total pressure, 
leaving the electron and ion temperature initially uni- 
form throughout the system. The density is given by 
([8]) and the magnetic field by ([9]), using the shape factor 
(p!Q|) , centered at ?/o, the middle of the domain in the y 
direction, with a half- width A = 0.5. A guide field of am- 
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FIG. 3. Ex electric field for run Snd (top panel), run Sr^ 
(middle panel) and run Su (bottom panel). The snapshot are 
averaged over around the same time as in fig [H 




0.02 



0.00 



0.02 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
10 20 30 40 50 60 

FIG. 5. Time evolution of the reconnect ion rates for the runs 
(black solid), S^d (green), Sr^ (dashed red) and Su (solid 
blue). 
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FIG. 4. Ey electric field for run S^d (top panel), run Sr] 
(middle panel) and run Su (bottom panel). The snapshot are 
averaged over around the same time as in fig[Tl 



plitude Bgf = 1 has been added for more generality. Al- 
though we do not show them here to be more concise, we 
have performed other test simulations within a coplanar 
initial condition and have observed a very similar behav- 
ior. The initial plasma temperature is T = 3/2 and the 
ion to electron temperature ratio is Ti/T^ = 5. Particles 



are loaded as locally Maxwellian distribution functions. 
This initial condition is again not a Vlasov equilibrium 
and the current sheet initially evolves slightly toward a 
state closer to a self-consistent equilibrium while emitting 
some ion scale waves in the system. These waves have 
not been found to affect the results discussed here in any 
way, although a kinetic steady state would certainly be 
preferable for both physical and practical reason^. This 
initial phase of asymmetric reconnection modeling will be 
addressed in further detail in a forthcoming study. 

n{y) = l-l[s{y) + S{yf) (8) 
B4y) = l + S{y) (9) 
5(2/) = tanh(^^^) (10) 

As in the previous section, we performed simulations 
with grid scale reconnection, resistive dissipation and 
hyper-resistive dissipation and investigate whether this 
changes lead to substantial differences in the dynamical 
behavior of reconnection. Runs with only numerical dis- 
sipation will be labelled A^'^, runs with classical resistive 
dissipation will be labelled A^'^, and those with hyper- 
resistive dissipation will be labelled A^'^. For each case, 
we perform two simulations with two mesh resolutions. 
The following analysis is performed at times occurring 
before t = r^, where ta = Xm/VA is the shortest Alfven 
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TABLE II. Summary of the parameters for the asymmetric 
runs 



Name 


r\ 


V 


l\y 






2 10"^ 





0.15 


300. 


^nd 


2 10"^ 





0.076 


150 


Aa 


1.4 10-2 





0.15 


150 




1.4 10-2 





0.076 


150 


Al 


2 10-^ 


5 10-^ 


0.076 


150 


Al 


2 10-^ 


5 10-^ 


0.05 


100 



travel time across the system. This guarantees that the 
periodicity of the system does not alter significantly our 
results. Table HIl summarizes the run parameters. 

A. Reconnection with numerical dissipation 

In this section, we analyze the results of runs A^'^. In 
these runs, the resistivity is so small that the reconnec- 
tion of field line is enabled by numerical dissipation. The 
two runs differ by the mesh resolution. The resolution, 
in the y direction, used for the runs A^^ and are 
At/ = 0.15 and A?/ = 0.076, respectively. In the x di- 
rection, both runs have the same resolution Ax = 0.15. 
These values are a factor of \/3 smaller in terms of the 
ion inertial length based on the tenuous side of the cur- 
rent sheet. Consequently, the mesh is sufficiently fine to 
resolve the scales at which the Hall effect becomes impor- 
tant. To save computer time, the downstream length of 
the domain used for run A^^^ is set to 150 while it is 300 
for run A^^. The Alfven time is about t\ ~ 115 for run 
and t\ ^ 60 for run A^^. The time step is At = 10"^ 
and At = 5 10"^ for A^^ and A^^, respectively. Such 
small time steps are necessary to prevent the violation 
of the CFL condition by the propagation of short wave- 
length whistler waves in the system. The width of the 
domain in the y direction is set to ym = 50 for both runs. 
There are approximately 264 millions macroparticles for 
each simulations. 

Figure [6] shows the out-of-plane electron bulk flow in 
simulation A^^ at times t = 50 (top panel) and t = 90 
(bottom panel). Surprisingly, the system starts with a 
very long period of time over which the current sheet 
elongates, while the flux tubes already reconnected by 
the initial perturbation are moving downstream. Then 
the reconnection process starts a globally unsteady phase 
and produces numerous plasmoids. After being born in 
the electron current sheet, these plasmoids grow quickly 
to ion scales and are slowly convected downstream. They 
dominantly grow on the weak field side owing to the 
weaker magnetic tension force there. Some plasmoids are 
ejected faster than others and coalesce with a preceding 
island. 

As mentioned previously, plasmoids are seen in fully 
kinetic simulations of symmetric systems and are inter- 
preted as the consequence of secondary tearing in long 
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FIG. 6. Out-of-plane electron bulk flow Vez for the run A^^ 
at t = 50 (top) and t — 90 (bottom). The in-plane magnetic 
field lines are represented by the solid white lines. 



and thin electron current layers. Although it has dif- 
ferent physical origins, this behavior is also reported in 
high-Lundquist number MHD simulations^^. To under- 
stand further the origin of the plasmoids in our simula- 
tion, we therefore focus on the length of the electron cur- 
rent layer. To actually measure this length, we take the 
sum of the distances between the dominant X point and 
the points where the separatrices are distant from each 
other by more than a given threshold, here chosen equal 
to 0.7 Si. Note that this measurement is not intended to 
be a quantitative measurement of the length of the cur- 
rent layer in itself, since it depends both on the arbitrary 
threshold and the angular aperture of the separatrices. 
It is however a simple measurement and a good quali- 
tative proxy of the elongation of the current sheet and 
the resulting length visually corresponds to what would 
have been identified as the edges of the layer. The black 
curve on Fig. [Tlis the result of this measure. As one can 
clearly see, the length of the layer increases rapidly to 
reach very high values and suddenly decreases at t ^ 50 
after the formation of multiple plasmoids (also seen on 
the top panel of Fig. [6]). It is interesting to notice that, 
all along the simulation, the current sheet does not stay 
short unless strongly influenced by two surrounding plas- 
moids. Their fast growth rate enables them to keep the 
current sheet short for a while before they move down- 
stream. As they move away, the current sheet starts 
to elongate again with roughly the same speed until in 
breaks once again. Although the data shown here stops 
at t = 100 because of a possible influence of the period- 
icity after the first Alfven time, the investigation of later 
times revealed similar behavior. The overall dynamics of 
magnetic reconnection in that system is thus dominated 
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by an oscillation between an elongation of the electron 
current sheet and its breaking in multiple islands. 

The blue curve on Fig. [71 represents the length of the 
electron current layer for the run A^^, where the mesh 
size has been divided by a factor 2. We observe the same 
behavior although this time the current sheet elongates 
much less than in the previous case. Assuming the plas- 
moids are the result of a secondary (and numerical) tear- 
ing instability, the observation of shorter current sheets 
can be understood as a consequence of the use of a higher 
resolution mesh. Because the dissipation is explicitly ne- 
glected in the equations, the current sheet collapses to 
the first and only non-ideal scale, the grid scale. In the 
high resolution run, the current sheet will therefore be 
thinner than in the coarser case. Accordingly, the cur- 
rent sheet does not need to be as long as in the coarse 
case to reach an aspect ratio where a numerical tearing 
is triggered and starts producing plasmoids. As a con- 
sequence of these shorter layers, instabilities occur more 
frequently but produce less plasmoids each time. One 
can conjecture that a run with a much higher resolu- 
tion could then trigger a numerical tearing instability for 
very short electron layers. Depending on whether these 
islands manage to reach large scales or are quickly recon- 
nected to exhaust field lines as they move downstream, 
the overall process could appear highly unsteady or spu- 
riously steady macroscopically, with a very short electron 
layer. 




FIG. 7. Length of the electron current sheet as a function of 
time for the run (black curve) and run A^^ (blue curve). 



rate of the higher resolution run is much larger than 
the one of the coarse run, but also much less steady. 
In particular, strong oscillations can be observed around 
t ^ 40, t ~ 50, t ~ 55. These times corresponds to 
times where strong variations of the current sheet length 
are also observed. This suggests that the frequent for- 
mation of plasmoids, by maintaining the current layer 
short, plays an important role in controlling the recon- 
nect ion rate. Whenever the formed plasmoids move far 
enough from their original location, the current sheet is 
less constrained and starts to elongate, and a decrease the 
reconnection rate is simultaneously observed (roughly, in- 
tervals t e [30,40], t e [43,48], t e [50,55]). The very 
fast triggering of a plasmoid instability restores a faster 
rate, as opposed to the coarser resolution run. In the 
coarser case, the current layer does not break up into suf- 
ficiently small current sheets to maintain a fast rate. It 
is interesting to notice that this overall behavior is phe- 
nomenologically quite similar to interpretation recently 
proposed in symmetric fully kinetic model^, although 
here the instability must have a numerical origin. 
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FIG. 8. Top : Time evolution of the reconnection rates for 
the runs (orange), A^^ (blue), A^^ (black) and A^ (green). 
Bottom : Reconnection rate for all asymmetric runs, aver- 
aged between t = 30 and t = 40, as indicated by the dash- 
dotted vertical lines on the top panel. 



The top panel of Fig. [8] shows the time evolution of 
the reconnection rate for the runs A^^ (green) and A^^ 
(blue). The two curves are significantly different. The 



We have performed two other simulations with the 
same resolution as run A^. but with smaller domain 
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sizes in the downstream direction. In all cases, we clearly 
saw an elongation of the electron current sheet but the 
following evolution was interestingly quite different. In 
the smaller case, where Xm = 64, we saw no plasmoid, 
whereas in the case where Xm = 150, we saw the current 
sheet to break only once in multiple plasmoids and, af- 
ter t ~ 60, stay short and localized at a single X point. 
Because the downstream size of the domain is the only 
difference between these runs, these features are the con- 
sequence of the periodicity influence on the system. In 
the smallest domain, the recirculation of plasma jets and 
electromagnetic fields occurs very soon. By artificially 
thickening the current layer, it prevents it from being un- 
stable to a secondary tearing and plasmoids thus never 
form. In the intermediate size system, the same effect oc- 
curs although the system size being larger, it allowed for 
the first set of plasmoids to exist but does not let them 
enough time and space for the current sheet to have the 
opportunity to elongate again. The current sheet there- 
fore stays short and the exhaust open, as a consequence 
of the periodic boundary and not because of a local mech- 
anism. We observed the elongation/breaking mechanism 
to stop at t ^ 60, which approximately corresponds to 
one Alfven time for such system, consistently suggesting 
an influence of the domain periodicity. This influence 
might be less important for higher resolution runs since 
the current sheet is then always shorter. 



B. Reconnection with resistive dissipation 

As for the symmetric simulations, we now explore the 
consequences of increasing the value of the Joule resistiv- 
ity up to 7^ = 1.4 10~^. The size of the domain is set to 
150 in the downstream direction, and 50 in the upstream 
direction. The resolution is set to Ax = 0.15 in both 
runs, while Ay = 0.15 and Ay = 0.076 for runs and 
A^, respectively. The time step is set to At = 10~^ for 
both simulations. There are approximately 132 millions 
particles in each case. 

In the case of symmetric reconnection, increasing the 
resistivity to this value had no big overall effect on recon- 
nection. The observed decrease of the reconnection rate 
was associated with the enhancement of large scale dif- 
fusion of electromagnetic energy and a small lengthening 
of the electron current layer. Surprisingly, increasing the 
resistivity in the asymmetric case now have important 
consequences regarding the large scale evolution of the 
system. As can be seen on the top panel of Fig. [9l made 



from run A^, the electron current sheet is short and lo- 
calized. A visual inspection of the layer at different times 
throughout the simulation (not shown) revealed that this 
snapshot is a good representation of what occurs in the 
system at other (later) times. The run gave almost 
exactly the same result and is thus not shown here. As 
can be seen on the left panel of Fig. [TOl the current sheet 
keeps approximately the same thickness when the resolu- 
tion is doubled, which indicates that the Joule diffusion 




FIG. 9. On all panels, the color represents the out-of-plane 
electron bulk velocity Vez and the white solid lines are the 
in-plane magnetic field lines. Top : run A^ at time t — 
50. Middle : Run A^ at time t = 34 when the amount of 
reconnected flux is identical to the top panel. Bottom : Run 
At at time t = 50. 
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FIG. 10. Jz measured along the y direction through the X 
point and averaged between t = 30 and t = 40 to reduce nu- 
merical noise. Left : Runs (blue) and A^ (black). Right 
: Runs A^ (blue) and (right). 



dominates the overall dissipation at the X line. However, 
the higher resolution case has a slightly higher current 
density, indicating again the difficulty of Joule diffusion 
to dissipate the incoming magnetic flux. The black curve 
on the top panel of Fig. [8] represents the reconnection 
rate for the run A^. Consistently with the steady aspect 
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of the current sheet, the reconnection rate is found to be 
rather constant, with a value around ~ 0.025. On the 
bottom panel of the same figure, one can see that both 
and have roughly the same reconnection rate. 



C. Reconnection with hyper-resistive dissipation 

Our last test consists in replacing the resistive term 
by a hyper-resistive term. In the symmetric case, hyper- 
resistivity was found to slightly decrease the rate com- 
pared to numerical dissipation, but the large scale fea- 
tures were observed to be unchanged, as opposed to the 
use of uniform Joule dissipation. In this asymmetric case, 
changing the dissipation physics again results in a dif- 
ferent macroscopic behavior. The right panel of fig. [3 
shows that for this value of hyper-resistivity u = 5 10~^, 
changing the resolution by a factor of 2 leaves the current 
sheet unaffected, as opposed to the case of Joule diffusion 
where a small different can be noticed. The structure of 
the current layer is also significantly different and con- 
trary to the resistive case, does not consist in a single 
peak anymore. Figure [9] shows the out-of-plane electron 
bulk velocity Vez for times t = 34.5 (middle panel) and 
t = 50 (bottom panel). The middle panel corresponds 
to a reconnected fiux equivalent to the one of the top 
panel, in the resistive case A^. At this same phase of the 
process, the two snapshots look globally similar, the elec- 
tron current being locally enhanced at the reconnection 
site in both cases and the field lines dominantly bent on 
the weak field side. A more careful examination reveals 
differences. At the reconnection site as well as in the sep- 
aratrix regions, the electron bulk fiows are weaker in run 
Arj than they are in run A^, . The left-right asymmetry re- 
sulting from the presence of an initial guide field seems to 
be more pronounced in the hyper-resistive run. Finally, 
the bottom panel of Fig. [9] shows that, at the same time, 
the reconnection process is much more advanced in the 
hyper-resistive case than it is in the resistive one. Con- 
sistently, the reconnection rate for the run A^'^ is much 
larger than for A^'^, as can be observed by comparing 
the black and orange curves on the top panel of fig. [9l 
The bottom panel shows that the higher resolution run 
A^ has the same reconnection rate. 



V. SUMMARY AND DISCUSSION 

We have performed two dimensional hybrid simula- 
tions to study the role of the dissipation scale physics 
in the process of magnetic reconnection and extend it 
to the more general case of asymmetric current layers. 
Standard hybrid models include ion kinetic physics and 
fiuid electron physics but neglect electron inertial effects 
(thermal and bulk). These effects are indeed small 
at the ion scales, and are usually thought to have no 
dynamical role beyond violating the frozen-in condition. 
With the hybrid model, the use of different non-ideal 



electron mechanisms allowed us to investigate whether 
the nature of the dissipation occurring at the X line can 
significantly alter the large scale reconnection dynamics 
or not. In this paper, we have used numerical, resistive 
and hyper-resistive dissipation. 

The first part of this study focused on a symmetric 
configuration, for which we have provided evidences 
that our model recovers basic features of 2D collision- 
less reconnection. The structure of the magnetic and 
electric Hall fields was found to be very similar to 
those measured in previous hybrid and fully kinetic 
calculations. Neglecting all dissipation mechanisms in 
Ohm's law, we observed the current layer to collapse 
down to the grid scale in both upstream and downstream 
directions, independently of the mesh size. However, as 
observed in previous studies, the overall reconnection 
rate seems unaffected and the process stays fast and 
steady, supporting the Hall reconnection paradigm. 
The use of a uniform resistivity increases the length 
of the current layer and also dissipates large scale 
electromagnetic energy in the exhaust region. As a 
result, the reconnection rate is significantly decreased. 
Peak amplitudes of electromagnetic fields at the ion 
scale are however unaffected by the viscous dissipation 
caused by hyper-resistivity, which remains at small 
scale, as previously understood. The reconnection rate 
is found to decrease a bit in comparison to the dissi- 
pationless case, however the difference remains small. 
A major difference between these results and modern 
fully kinetic calculations is the lack of plasmoids. In all 
cases here, the reconnection process was indeed steady. 
Whether this discrepancy means that hybrid models 
lack key physical ingredients is however still unclear, 
the formation mechanism of these plasmoids and its 
dependance on the inherent and over-estimated shot 
noise in Particle-In-Cell simulations being still unknown. 
Furthermore, if in situ measurements have provided evi- 
dences for the observation of such magnetic island o^^i^^ , 
the statistics of their occurrence his however unknown, 
leaving unanswered the question of their importance 
regarding the overall steadiness of magnetic reconnection 
in collisionless systems. Fully kinetic Vlasov simulations 
can address this issue but remain computationally 
challenging considering the large domains required. 

In a second part, we have implemented an asymmet- 
ric initial condition, already used in previous works and 
again performed several calculations with different non- 
ideal mechanisms. Contrary to the symmetric case, the 
nature of the non-ideal mechanism violating the frozen- 
in condition has been found to lead to substantial differ- 
ences regarding the macroscopic reconnection dynamics. 
When the dissipation is neglected, the electron current 
layer is seen to elongate until it triggers a plasmoid in- 
stability. Because they grow faster than they move, the 
plasmoids constrain the current layer, which stays short 
for a while but starts elongating again as the islands move 
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downstream. Consistently, a higher resolution simula- 
tion triggers the instability for shorter layers and leads to 
more frequent instabilities producing less plasmoids. As 
a result the current layer stays dynamically shorter and 
the reconnection rate is globally faster, but very unsteady 
because modulated by the island formation process. 

Surprisingly, the use of a uniform and quite large re- 
sistivity was observed to prevent the electron layer from 
elongating and no plasmoids were observed in such cases. 
Similarly, the hyper-resistivity keeps the electron layer 
short and the process steady. Like in the symmetric 
case, controlling the thickness of the current layer with a 
uniform resistivity without enabling large scale diffusion 
appears difficult. 

This work suggests two major questions : 1/ Why 
is the electron layer elongating in the absence of phys- 
ical dissipation in the Ohm's law while the Hall term is 
present? 2/ Why does it occur only in the asymmet- 
ric case? Let us consider a case where fast reconnection 
occurs in a steady state, and at a given time the dissi- 
pation mechanism ceases for some reason. Accordingly, 
the field lines, in the upstream region, that were about to 
be reconnected cannot be reconnected anymore. In the 
exhaust however, the field lines that have already been 
reconnected will keep moving downstream and carry with 
them the hot plasma away from the reconnection region. 
As a result of these two effects, a pressure unbalance will 
force the localized electron layer to elongate and become 
a ID tangential layer. If one turns the dissipation mech- 
anism on again, new reconnected flux is provided to the 
exhaust, along which Hall disturbances can propagate 
and open the exhaust again, increasing the reconnection 
rate. This interpretation is similar to the localizing effect 
of the non-ideal electron term, recently proposedi^i^. If 
numerical dissipation can violate the frozen-in condition 
and enable reconnection of field lines, one cannot, how- 
ever, expect that it behaves like an actual dissipative 
effect, i.e. there is no reason for which it should have 
this localizing property and therefore an elongation of 
the current layer can be expected. 

In the symmetric system however, the electron current 
sheet stays localized even in the absence of physical 
dissipation in the Ohm's law, which seems to contradict 
our previous interpretation. Let us notice the following 
facts: in this system, the exhaust is widely open on 
both upstream sides of the neutral line, as opposed 
to the asymmetric case where the field lines on the 
strong field side are barely bent. In the asymmetric 
dissipationless runs, we have noticed that two plasmoids 
growing rapidely around the X point can, for some time, 
maintain the electron layer short. Finally, we have also 
observed that the widening of the exhaust region due 
to the periodicity of the domain and the associated 
recirculation of plasma, can stop the elongation of the 
current layer and keep a single, localized, X point. 
Although further work is needed to fully understand this 
effect, these facts already suggest that the dissipation 
occurring at the X line may not be the only mechanism 



acting to constrain the length of the electron layer, and 
a feedback from the ion scale, whether it is controlled 
by the Hall effect, the growth of a nearby plasmoid or 
the widening effect associated to periodic recirculation, 
might also play an important role. The asymmetry of 
the initial system strongly affects the ion scale geometry 
surrounding the X line, which might then lead into a 
different constraint on the length of the electron current 
sheet than in symmetric systems. 



In all the simulations presented in this work, the field 
line connectivity can be changed, whether this is due 
to grid effects, resistivity or hyper-resistivity. However, 
not all simulations show with the same overall evolution. 
Therefore, this work provides evidences that collisionless 
reconnection is generally not independent from the dis- 
sipation effects, as if they could be anything as long as 
they enable the change of connectivity. It is however pos- 
sible that, for some reason that is yet not understood, the 
coupling of the Hall term with the appropriate dissipa- 
tion term results in an evolution that, apparently, does 
not depend, or very weakly depends, on the parameter 
controlling the dissipation. Further work is needed to un- 
derstand to what extent the reconnection rate depends 
on this localized dissipation once it is operating. This is 
beyond the scope of this paper and will be the topic of 
a future study. Preliminary results suggest that varying 
slightly the hyper-resistivity coefficient does not change 
the reconnection rate. Considering that, for these tests, 
smaller values of hyper-resistivity would be preferred to 
larger ones, the interval of values that can actually be 
tested is limited by the high resolution and small time 
steps required. These, with the requirement of long do- 
mains in the downstream direction, make such study dif- 
ficult. Although hyper-resistiviy shares common proper- 
ties with the kinetic dissipation in symmetric systems, 
it remains a simplified model and there is no theoreti- 
cal proof that the similarity persists in asymmetric sys- 
tems. Furthermore, the link between the non-uniform 
and self-consistent kinetic dissipation coefficient obtained 
from collisionless theory and the value of the uniform 
hyper-resistivity is not clear even in symmetric systems, 
therefore comparing the variation of the hyper-resistivity 
with the variation of the electron mass in a fully kinetic 
system is not trivial. Understanding the kinetic mech- 
anisms leading to dissipation in asymmetric systems is 
thus an important issue and is crucial for the under- 
standing of the future NASA Magnetospheric MultiScale 
(MMS) mission. Preliminary results in the comparison 
between hyper-resistive hybrid and fully kinetic PIC sim- 
ulations of the same asymmetric systems reveal striking 
similarities and a detailed comparison is the topic of a 
separate paper—. 
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